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Abstract 

We present a method to perform a least squares fit of a decay chain involving mul- 
tiple decay vertices. Our technique allows for the simultaneous extraction of decay 
time, position and momentum parameters and their uncertainties and correlations 
for all particles in a decay chain. 



1 Introduction 



In high energy physics experiments decay reactions that proceed via inter- 
mediate metastable states are usually reconstructed by following a bottom-up 
approach. One starts by extracting the parameters of those decay vertices from 
which the reconstructed final state particles emerge and uses the intermedi- 
ate 'composite' particles for the reconstruction of upstream decays. At each 
decay vertex the parameters of the composite particle are determined with a 
least squares fit to its daughter particles, subject to the constraint that those 
originate from a common point. The disadvantage of this approach, which is 
sometimes called 'leaf-by-leaf ' fitting, is that constraints that are upstream of 
a decay vertex do not contribute to the knowledge of the parameters of the 
vertex. An example of a decay for which this is impractical is — > tt^tt". 

In this paper we discuss the implementation of a least squares fit that extracts 
all parameters in a decay chain simultaneously. We shall call this fit, which 
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wc developed for data analysis in the BABAR experiment, a global decay chain 
fit. First, we propose a parameterization of a decay chain in terms of vertex 
positions, momenta and decay times. Subsequently, we argue that the Kalman 
filter is a suitable technique to extract these parameters and the corresponding 
covariance matrix from the external constraints, which in the case of BABAR 
are reconstructed charged particle trajectories and neutral particle calorimeter 
clusters. Finally, we present two examples and briefiy summarize experience 
with the fit in BABAR. 

The decay chain fits discussed here are hypothesis driven. The task of finding 

the reconstructed tracks and clusters and associating those with the final state 
particles in the decay tree is outside the scope of this paper. In BABAR physics 
analyses decay trees are built layer-by-layer, usually by making all possible 
combinations of final state particles and applying selections on the invariant 
mass and vertex x^. Vertex pattern recognition plays an insignificant role 
because the low combinatorics does not warrant more complicated algorithms 
and because the track parameter resolution is barely sufficient to separate the 
decay vertices of the particles that are of most interest to the experiment, 
namely B and D mesons. 



2 Pcirameterization of a decay tree 

Figure 1 shows a schematic picture of a decay tree. The positions of the ver- 
tices in the decay tree, and the momenta of all particles, constitute the degrees 
of freedom of the decay tree. These degrees of freedom, the internal constraints 
(such as momentum conservation at each vertex) and the relation to the ex- 
ternal reconstruction objects, define the decay tree model. 

The choice of parameters in the decay tree model is not unique, but we found 
the following parameterization suitable for use in BABAR. Each reconstructed 
or 'final state' particle is represented by a momentum vector {Px,Py,Pz)- The 
mass of a final state particle is not a parameter in the fit, but assigned based 
on the particle hypothesis in the decay tree. Each intermediate particle in 
the decay tree is modeled by a four momentum vector {px,Py,Pz, E) and a 
decay vertex position (x,y,z). If the composite particle is not at the head of 
the decay tree, we also assign a parameter for its decay time. We choose this 
parameter to be ^ = where / is the decay length. 

If a composite particle has an expected decay length much smaller than the 

vertex detector resolution, we call this particle a 'resonance'. A resonance 
does not have a decay time parameter and it shares the decay vertex position 
with its mother, unless it is at the head of the decay tree. In the BABAR 
reconstruction software particles with an expected decay length cr < 1 /im, 
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Fig. 1. Schematic picture of a decay tree with three charged particles reconstructed 
as track segments (T), one photon reconstructed as a calorimeter cluster (C), and 
two composite particles (I for 'internal' and H for "head'). 



such as 7r°, J/V^ and D* ^ are treated as resonances. ^ 

We distinguish two types of constraints in the decay tree. Two internal con- 
straints are applied to remove redundant degrees of freedom: the vertex con- 
straint expresses the relation between the decay vertex of a particle and the 
production vertex of its daughters; the momentum conservation constraint 
ensures four-momentum conservation at each vertex. The reconstructed final 
state particles constitute the external constraints. In this paper we consider 
only 5-parameter track segments and calorimeter clusters with a reconstructed 
position and energy. Exphcit expressions for the constraints are given in sec- 
tion 4. To put those in the proper context we introduce the fit procedure 
first. 



We use the term resonance for any particle with a short lifetime, regardless of 
whether the decay is through the strong or electroweak interaction. 
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3 Fitting a decay tree 



3.1 Measurement constraints and exact constraints 

To extract the optimal set of decay tree parameters from the reconstruction 
information we use a least squares fit. Within this framework constraints are 
expressed as contributions which are a function of the parameters of the 
model, collectively denoted by the vector x. The solution to the fit is the 
value of X that minimizes the total x^- To define the contributions we 
distinguish between exact constraints and measurement constraints. The latter 
are characterized by the fact that they have an associated uncertainty, whereas 
the former do not. The internal constraints of the decay tree model are exact 
constraints, whereas the external constraints are measurement constraints. 

For a measurement constraint i the contribution takes the form 



where r.;(x) is the constraint residual and Vi is the constraint variance. We use 
a matrix notation so that the x^ contribution is defined for constraints of any 
dimension. The residual can be defined as 



where is the value of the measured quantity (for example the parameters 
from a reconstructed track) and hi{x) is the measurement model that expresses 
this quantity in terms of the parameters x. With this definition Vi is the 
variance in the measurement m^. 

Exact constraints can be imposed on the model by parameter substitution. For 
example, the parameters that represent the momentum of a mother particle 
can directly be replaced by a sum over the daughter momenta. Although this 
simplifies the fitting problem by reducing the number of parameters, it is not 
suitable for the decay tree fit discussed here. First, the substitution procedure 
makes it more complicated to use recurrence in the implementation of the 
fitting algorithm. Second, it complicates the calculation of those parameters 
that have been removed — the momentum of the mother — , especially where 
it concerns the associated uncertainty. 

Consequently, exact constraints are not applied by substitution but by means 
of Lagrange multipliers. If the exact constraint in terms of x is written as 



xl = rj{x) V; 1 ri{x) 



(1) 



ri{x) 



rrii - hi{x) 



(2) 



9i{x) = 0, 



(3) 
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the contribution of the constraint takes the form 



= 2Xjg,{x) 



(4) 



where Aj is the Lagrange multipher. The latter is added to the hst of param- 
eters X for which we minimize the x^. 

3.2 The standard least squares fit 

The requirement that x^ = J^iXi be minimal defines a set of equations 
dx^/dx = 0. A solution can be obtained with a Newton- Raphson method. 
(See for example [1].) Given an initial estimate x^^\ one obtains a new esti- 
mate 



For constraints that are not linear in x this expression can be applied iter- 
atively until a certain convergence criterion is met. Using error propagation 
one can derive that the variance in x is given by 



Wc shall call the procedure described above, in which the x^ contributions of 
all constraints are minimized simultaneously, the standard fit. ^ The expres- 
sions for the standard fit show that the minimization procedure requires the 
inversion of matrices with the dimension of the parameter vector x. Compli- 
cated decay chains can easily require several tens of parameters, which leads 
to large computational costs. 

3.3 The progressive least squares fit 

The Kalman filter or progressive fit [2] is a x^ minimization method that is 
less computation intensive than the standard fit. In the reconstruction of data 
from particle physics experiments the Kalman filter is mainly applied in track 
fitting, where its virtue is not only speed, but also the possibility to easily 

^ This method for x'^ minimization is often called a 'global' fit, but we reserve the 
term global in this paper for use in 'global decay chain fit', a decay chain fit that 
determine all parameters and correlations in a decay tree simultaneously, rather 
than leaf-by-leaf. 




(5) 




(6) 
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include multiple scattering effects as so called 'process noise' [3]. The Kalman 
filter has been proposed for vertex fitting by several other authors [3.4]. How- 
ever, a general treatment of exact constraints has not been considered before, 
which justifies a re-derivation of the Kalman filter for the purpose of decay 
chain fitting. To connect with previous work we follow a notation close to that 
in [3]. 

Consider a measurement constraint k. Assume that the consisting of all 
contributions of constraints {0, . . . , A; — 1} has already been minimized with 
respect to the model parameters, leading to a 'prediction' Xk-i with variance 
Ck-i- To calculate how x changes when the constraint k is included in the 
minimization, we define a contribution 

xl = (x - Xfc-i)^ Cfc_r^ (x-Xk-i) 

+ {hk{x) - rukf Vk~^ {hk{x) - nik). (7) 

The least squares solution for x is obtained from the requirement that dx^/dx = 
0, i.e. 

Ck-i-\x - Xk-i) + Hi {hk{x) -nik) =0 (8) 

where H^. = dh/dx\x^^_^ is the derivative or projection matrix. We call the 
solution to this equation the updated parameter vector x^. 

The authors of reference [3] discusses two approaches to calculate Xk, which 
are called the 'gain matrix formalism' and the 'weighted means formalism'. 
The latter still requires inversion of matrices with the dimension of x and is 
not suitable for our purposes. In the gain formalism one rewrites Eq. 8 as 

(Cfe_r' + HlVk-'Hk) {x - Xk-i) = HlVk-'rt' (9) 
where we defined the so-called residual of the prediction 

r-fe"^ = mfe - hk{xk-i) (10) 
and where we assumed that the measurement model is linear, i.e. 

hk{x) = hk{xk-i) + Hk{x - Xk-i). (11) 

The treatment of more general h will be discussed below. Solving for x, we 
obtain 

Xk = Xk-i + Kk rl-^ (12) 
where the gain matrix Kk is defined as 

Kk = (Ck-i-^ + HlVk-^Hk)'^ HlVk-\ (13) 
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The latter can be rewritten as 

K, = Ck-,Hl(Rl-'y' (14) 
where 

^ Vk + HkCk-iHl (15) 

is the uncertainty in the predicted residual. The fact that Eq. 14 contains 
the inverse of a matrix with the dimension of the measurement m^, rather 
than with the dimension of the parameter vector ,t, is the reason that the 
progressive fit is in general faster than the standard fit. 

The updated covariance matrix can be obtained by error propagation from 
Eq. 12, which gives 

Cfe = (1 - KkHk) Ck-i (1 - KkHkf + KkVuKu". (16) 



This expression is computation intensive because the first term on the right 
hand side corresponds to a product of three square matrices with the dimension 
of the parameter vector. It can be simplified to 



Cfe = (1 - KkHk)Ck-i (17) 

which is much faster, but known to be sensitive to finite machine digit ac- 
curacy, in particular if Vk is small compared to HkCk-iH^ [5]. This can be 
understood by evaluating the effect of a small perturbation in the gain ma- 
trix Kk ^ Kk + SK. Substituting this in Eq. 16 yields 5Ck — SKR^SK^ , 
whereas Eq. 17 gives 5Ck = —SKHkCk- As a result the second expression can 
lead to a covariance matrix with a negative determinant, an effect that we in- 
deed observed in fits with many parameters. We have found that by rewriting 
Eq. 16 as 

Cfe = Ck-i — Kk (2HkCk-i — RI~^K]^^ , (18) 

the computational stability is preserved at the cost of a moderate increase in 
computation time with respect to Eq. 17. 

Finally, the contribution of the constraint ruk is given by the contribution 
of the prediction, i. e. 

xi - irtriRtr'rr- m 

The contribution of a particular constraint depends on the order in which 
the constraints are applied. However, if all constraints are linear in x, the sum 
of the contributions is independent of that order. 
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3-4 Non-linear constraints 



If the measurement model is not linear in x, Xk can be extracted with an 
iterative procedure. Expanding h{x) around a point x^*^ (initially given by 

Xk-l), 

hk{x) = + H'i\x-x^'^) (20) 

the estimate for Xk becomes 

xt^'^ = X,., + rr^« (21) 
where the residual ^^^^ is given by 

rJ-^W = m, - /i,(x«) - Hi'\x,_^ - x«). (22) 

Choosing the new expansion point, one obtains an improved esti- 

mate of Xk by iteration, subject to a certain convergence criterion. A suitable 
observable to test the convergence is the contribution 

The time consuming calculation of the covariance matrix (Eq. 16) can be 
performed after convergence is obtained. 



3.5 Exact constraints 



The contribution for an exact constraint g{x) = was introduced in Eq. 4. 
Minimization of the under the exact constraint is performed by solving the 
set of equations dx^/dx — and dx^/dX = for x and A simultaneously. To 
derive the expressions for the progressive fit we define a x^ 

Xl = {x - Xk-ifCk-i-\x - Xk-i) + 2Xlgk{x). (24) 
Linearizing the constraint equation around x^-i 

9k{x) = gk{xk-i) + Gk{x- Xk-i) + ■■■ (25) 
where G — dg/dx, we obtain for the linearized set of first derivatives 

= Ck-i-\x - Xk-i) + GlXk 

(26) 

= gk{xk~i) + Gk{x- Xk-i) 
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Multiplying the first equation by GCk-i and subtracting the second equation 
yields 

(GfeCfc.iCD Afe = gk{xk-i). (27) 

Since for non-trivial constraint equations the matrix on the left side is invert- 
ible, this equation leads to a solution for A^. Eliminating we obtain for the 
updated parameter vector 

Xk = Xk-i - Kkgk{xk-i) (28) 
where we defined the gain matrix for exact constraints by 

Kk — Ck-iGj. (CkCk-iGl 

T ■ (29) 
Using error propagation we obtain for the updated covariance matrix 

Cfe = (1 - KkGk) Ck-i (1 - KkCkf (30) 
and for the contribution 

Xk = gk{xk-i)'^ [GkCk-iGl^ ^ gk{xk-i). (31) 

A comparison of the equations for the exact constraint with those for the 
measurement constraint shows that with the substitutions 

h{xk-i) -ruk — > gk{xk-i) 

{o2) 

Vk + HkCk-\Hj^ — > GkCk-iG^. 

the two procedures are identical. This non-trivial result confirms the intuitive 
notion that an exact constraint is effectively the same as a measurement with 
infinite precision. 

The minus sign in equation Eq. 28 with respect to its counterpart Eq. 12 is 
the result of our decision to choose a conventional notation. It is customary to 
express the Kalman filter equations in terms of the derivate of h{x) to x rather 
than the derivate of the residual r{x) to x. The minus sign reflects the fact 
that these derivatives differ by a sign when the residual is defined as in Eq. 2. 
For the implementation of the fit we have defined the residual with opposite 
sign, such that the same expressions can be used for exact constraints and 
measurement constraints. 

It is interesting that the progressive fit deals more effectively with exact con- 
straints than the standard fit does. In the latter the Lagrange multipliers are 
added explicitly as parameters to the fit, increasing the number of parameters 
and therefore computational costs. The treatment of an exact constraint in 
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the standard fit is therefore relatively expensive. In the progressive fit the La- 
grange multipliers can be eliminated and exact constraints are less expensive 
than measurement constraints. 



4 Explicit expressions for constraints in the decay tree 

In this section we provide explicit expressions for the internal and external 
constraints introduced in section 2. 



4-1 The internal constraints 

The internal constraints reduce the set of decay tree parameters to a set that 
is not overcomplete. In the absence of a magnetic field or for neutral particles, 
the momentum vector p is constant, leading to a time-evolution of the position 
x{t) = xo + pt/'-fm, where t is the time in the lab frame, m the rest mass 
and Xo the position at t = 0. For each composite particle i in the decay tree 
the vertex constraint expresses the relation between its decay vertex Xi, its 
momentum pj and the decay vertex of its mother xm, 

xm + 9iPi = (33) 

where 9 = l/\p\ = t/'~fm is the decay time parameter. The momentum con- 
straint for particle i can be expressed as 

Pi-Y.Pj = and Ei-Y,Ej = 0, (34) 

3 3 

where the sum runs over the momenta of all the daughters. 

For charged particles in a magnetic field the expressions for the internal con- 
straints are more complicated. In most practical applications one can neglect 
the effect of the magnetic field for the particles inside the decay tree, because 
the bending in a typical lifetime is very small. This is for example the case for 
B and D mesons in the 1.5-T magnetic field of the BABAR detector. It is not 
true for some charged baryons such as and . The full expressions for 
the internal constraints in a constant magnetic field have been implemented 
for BABAR but are outside the scope of this paper. 
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4-2 The reconstructed track constraint 



The magnetic field in the BABAR spectrometer is approximately homogeneous 
and aligned with the 2;-axis of the detector. Reconstructed charged particle 
trajectories are represented by a 5 parameter helix m^^ = {do, (po,uj, zo,taii X) 
and a corresponding covariance matrix. In terms of these parameters, the 
position X and momentum p of the particle along the helix trajectory are 
given by 



X = r sin(0) — (r + do) sm{(f)o) 
y = -r cos(0) + (r + do) cos(0o) 
z — Zo + /tan(A) 

Px = Pt COS(0) 

Py =ptsin((/)) 
Pz = Pt tan(A) 



(35) 



where / is the flight length in the transverse plane, measured from the point of 
the helix closest to the 2;-axis, (p = (po + r = 1/ui and pt = qa/u with q the 
charge of the particle in units of the positron charge and a[J/m\ = —e[C]Bz[T] 
a constant. The measurement model h{x) used in the decay tree fit follows from 
the inverse transformation and can be written as 





( 4 \ 


' {pm-Pt)/aq ^ 






aX{m2{pyo,PM) 


h = 


UJ 


aq/pt 




zo 


z - Ipz/pt 




^ tan A J 


^ Pz/Pt 1 



(36) 



with Pt = ^pI+P% Pxo 
(p — atan2(py,pa;) and I - 



PxO + PyO-i 



p^ + aqy, Pyo = Py - aqx, pto = 
(p — (po)pt/qa. The derivatives can be concisely 
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written as 
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4-3 The reconstructed cluster constraint 



In BABAR reconstructed calorimeter clusters are represented by a measured 
position and energy = (^cius- ^cius) and a corresponding covariance matrix 
Vcius- Given a mother decay vertex x and a momentum vector p, the measure- 
ment model is defined as 



h 



y + Opy 

z + Qp^ 

\IpI + pI + p\ 



(38) 



where 9 takes the role of the photon 'decay length'. In principle, 9 can be added 
to the photon parameter list and extracted from the fit. Since this parameter 
is not very interesting, it is preferable to eliminate it. This can be done by 
reducing the set of four constraint equations r{x) = m — h{x) = to three by 
redefining the residual, for example 



^ (Xclus - X)py - (yclus - y)Px ^ 



r'{x) = 



('f'clus -^^Pz (^clus ^)Px 

^cius - ^JpI +pI+pI ) 



(39) 



The choice of r' is not unique: This particular choice is not suitable if pj. is zero, 
because the two position constraints would no longer be independent. However, 
for every value of p a set of independent constraints can be constructed. The 
variance of the new constraint is simply given by 



(40) 



where P — dr' /dm is the derivate of the constraint to the measurement. 
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4-4 Other constraints 



The constraints discussed above constitute the minimal set of constraints nec- 
essary for fitting a decay tree with final state particles that are reconstructed 
either as charged tracks or as calorimeter clusters. We have considered and 
implemented several other constraints. Sometimes, position and momentum 
parameters can be improved by constraining the mass of composites in the 
decay tree to the known particle mass. Knowledge about the interaction point 
can be used to constrain the production vertex of the head of the decay tree. 
Information about the beam momenta can be used to constrain the four- vector 
of the head of the decay tree, e.g. in e+e" T{4S) B^B^ decays. In addi- 
tion to charged tracks and photons, we have used reconstructed particles 
which are detected in the BABAR calorimeter, but for which the deposited 
energy is not useful as an estimate of the magnitude of the momentum. Fi- 
nally, we have found that missing particles can be included in the decay tree, 
provided the tree is not kincmatically under-constrained. The expressions for 
these constraints are easy to derive so we will not include them in this paper. 



4.5 Ordering constraints 

A disadvantage of the progressive fit with respect to the standard fit is that 
the final result of the fit can be sensitive to the order in which constraints 
are applied. In fits with 'process noise' such as track fits, there is a natural 
ordering. However, in fits without process noise, the ordering of the constraints 
can be chosen freely. 

The order in which constraints that are linear in x are applied is irrelevant, be- 
cause the covariance matrix contains all essential information on the constraint 
derivatives. If, for example, a constraint k is processed such that gk{xk) = 0, 
then gk{xn) = for any n > k. This is not the case for non-linear constraints, 
since the covariance matrix does not contain information on the higher deriva- 
tives of the constraint equation. Consequently, the order in which constraints 
are applied becomes important: The most non-linear constraints should be 
applied last. 

One can consider applying correlated non-linear constraints simultaneously, 
rather than consecutively. For example, the vertex constraint can be applied 
separately in the three space coordinates x, y and z. However, the dependence 
on the decay time parameter 9 correlates these constraints. Since they are also 
slightly non-linear, it is preferable to treat them as a single three-dimensional 
constraint rather than three separate one- dimensional constraints. Ultimately, 
all constraints can be applied simultaneously, which effectively reduces the 
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progressive fit to a standard fit. The performance advantage of the progressive 
fit then practically disappears. ^ 

Based on these considerations we have chosen the following approach to order- 
ing and combining constraint. The external constraints (reconstructed tracks 
and cluster, constraints to the interaction point etc.) are treated first. Subse- 
quently, all four-momentum conservation constraints are applied, starting at 
the end of the decay tree. Finally, at each vertex the geometric constraints 
and an eventual mass are combined. These combined constraints are applied 
consecutively starting at the end of the decay tree. 



4-6 Fit initialization and convergence 

The progressive fit requires initialization of both the parameters xq and the 
covariance matrix Cq. Vertex positions are initialized with the average interac- 
tion point or with the point of closest approach of (a subset of) reconstructed 
track segments in the decay tree. The momenta of particles reconstructed as a 
track segments are initialized by evaluating the track parameters at the point 
of closest approach to the initial vertex positions. The momenta of photons 
are initialized by using the initial vertex positions as their origin. Particle mo- 
menta inside the decay tree are initialized by adding the initial four vectors 
of their daughters. Finally, the decay time parameters are initialized from the 
initial vertex positions and momenta or from the expected decay time. 

The covariance matrix must be initialized with uncertainties that are large 
enough that the initial parameters have a negligible weight in the final result 
of the fit, yet small enough that the measurement errors Vk in equation are not 
numerically negligible with respect to HkCkH^- We have chosen for a diagonal 
matrix with diagonal elements that are roughly a factor 1000 larger than the 
square of the typical resolution for the corresponding parameter. 

Even with the iteration of non-linear constraints described in section 3.4 the 
decay chain fit does not converge in a single processing of all constraints. 
Therefore, we repeat the fit procedure until the total is stable. At each 
step xq is initialized with the result of the previous step, while Co is reset to 
its original value. Fits typically converge in three iterations. 



^ The progressive fit still has the advantage that exact constraints do not lead to 
extra (Lagrange multiplier) parameters. 
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5 Experience in BABAR 



We have implemented the decay chain fit described in the previous sections 
in the BABAR analysis framework. The fit has been extensively tested and is 
used in a variety of physics analyses. It is a useful alternative for leaf-by-leaf 
vertex fits, in particular for the reconstruction of decay chains with vertices 
with insufficient downstream constraints, such as —>■ A^n^ and Kg —>■ tt'^tt^, 
consistent treatment of mass constraints at several levels in a decay tree, direct 
extraction of the B^B^ decay time difference for analysis of time-dependent CP 
violation in T(45') — > B'^B^ events and kinematic fits with missing particles. 
We briefiy discuss two examples. 

Figure 2a shows the cr distribution of candidates in a BABAR simulation of 
B° — > D'^TT" decays with K'n'^n'^. Only candidates that are matched 

to the Monte Carlo truth are shown. The average reconstructed lifetime is close 
to the average lifetime cT]j+ = 0.311 mm with which the events were generated. 
Figure 2b shows the cr pull distribution, which has an RMS consistent with 
one. 




c x^, pull 



Fig. 2. Reconstructed lifetime (a) and lifetime pull (b) for extracted from a 
global decay chain fit to simulated 5° — > D+tt" (-D+ K^tt'^tt'^) decays. 
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As a second example, we study the reconstruction of — >■ 7r°7r° candidates 
in simulated J/iJjK^ decays with tt" 77 and J/ip — > fi^fi^. The 

7r°7r° decay cannot be reconstructed with BABAR^s traditional leaf-by- 
leaf fit. A global decay tree fit can be performed, provided that the origin of the 
Kg is known and that mass constraints are applied to the 7r° candidates. In this 
particular topology the J/ip vertex provides the origin and the B^ — > J/'^K^ 
fit is thus sufficiently constrained. 

Figure 3a shows the invariant mass distribution of the K^^ ti^-k^ candidates 
extracted from a fit to the B^ — > J/i/^Kg decay tree. The central value and 
resolution are significantly better than for the 'raw' tt^tt^ invariant mass dis- 
tribution that was obtained by assuming that the 7r° decays originate directly 
from the J/tJj — > jj^fi^ vertex. Figure 3b shows the consistency of the 
5° J/ipK'^ decay tree fit. (The has two degrees of freedom.) It is not 
entirely consistent with a flat distribution, partially because the energy of the 
photons originating from the 7r° candidates is not always fully reconstructed 
in the BABAR calorimeter. 
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Fig. 3. Kg — > vr'^vr*^ invariant mass (a) extracted from a global decay chain fit to 
simulated —>■ J/tpKg decays and the consistency of the fit (b). The hashed 
distribution in figure (a) is the tt^tt^ invariant mass before the geometric fit is 
performed. 
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The performance of the fit has been compared to a traditional leaf-by-leaf 
vertex fit that has been used in BABAR for several years. A direct comparison 
of computational performance is not trivial since that depends on implemen- 
tation choices in addition to the algorithmic complexity of the problem. The 
complexity of the global decay chain fit grows roughly with the square of the 
number of vertices in the decay chain, whereas the leaf-by-leaf fits (which do 
not keep track of correlations) behave more linear. Despite this fact, com- 
putation time use has not been a limitation in practical applications of the 
fit. 
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